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We derive constraints on the time variation of the Higgs vacuum expectation value {<j)) through 
the effects on Big Bang nucleosynthesis (BBN) and the cosmic microwave background (CMB). In 
the former case, we include the (previously-neglected) effect of the change in the deuteron binding 
energy, which alters both the '*He and deuterium abundances significantly. We find that the current 
BBN limits on the relative change in {(p) are —(0.6 — 0.7) x 10~^ < A{(j>)/{(/>) < (1.5 — 2.0) x 10~^, 
where the exact limits depend on the model we choose for the dependence of the deuteron binding 
energy on (0) . The limits from the current CMB data are much weaker. 
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I. INTRODUCTION 

Physicists have long speculated that the fundamental 
constants might not, in fact, be constant, but might in¬ 
stead vary with time P|. Among the various possibili¬ 
ties the interaction coupling constants have received the 
greatest attention (for a recent review, see Ref. 0). 
A great deal of attention has been focused on the fine 
structure constant a, for which a time variation has been 
claimed in observations of quasar absorption lines Q- 
Cosmological limits on a time variation in a can be de¬ 
rived from both big-bang nucleosynthesis (BBNland the 
cosmic microwave background (CMB) [4 IMIrF71HHIi(T|. 
Comparatively less interest has been inspired by possible 
time variation in the strength of the weak interaction, as 
parametrized by the Fermi coupling Gp. 

However, as emphasized by Dixit and Sher m , the 
Fermi constant is not a fundamental constant in the same 
sense as the fine structure constant. The Fermi constant 
is given by 

— = ^ ( 1 ) 

Therefore, it is more appropriate to examine the varia¬ 
tion of the Higgs vacuum expectation value than that of 
the Fermi constant. In fact, the time variation of the vac¬ 
uum expectation value of a field seems somewhat more 
plausible than the time variation of a fundamental cou¬ 
pling constant. (For arguments in favor of a possible 
spatial variation of the Higgs vacuum expectation value, 
see Ref. [l^l. 

Limits on the time variation of the Higgs vacuum ex¬ 
pectation value have been derived from BBN imm. 


However, improved observational limits on the primor¬ 
dial element abundances now allow us to place stronger 
limits than those derived, for example in Ref. [T^ . Fur¬ 
thermore, these previous studies ignored the effect of the 
change in the deuteron binding energy, which we have 
incorporated into our calculations. 

The effects of a variation in ((p) on the CMB have been 
investigated previously in Ref. |l4l| , but this paper made 
no attempt to derive actual limits based on observational 
CMB data, which we will do here. 

In this paper, we derive constraints on the time vari¬ 
ation of {(p) from both BBN and CMB observations. In 
this regard, the present work is most similar in spirit to 
that of Avelino et al. |^, who did a comparable calcula¬ 
tion for a. As noted in Ref. ^3, the effect on the CMB 
of changing the Higgs vacuum expectation value is sim¬ 
ilar, but not quite identical, to the effect of changing a. 
The effects on BBN, on the other hand, are quite differ¬ 
ent. We assume that {(jf) has the same value at the era 
of primordial nucleosynthesis T ~ 10^° — 10® K as it has 
at recombination T ~ 10^ K, but our results can easily 
be generalized to cases where this is not so. 

In the next section, we discuss the effects of {(p) varia¬ 
tion on BBN. In Sec. HI, we present the effect on recom¬ 
bination and calculate the CMB temperature anisotropy. 
Finally, we discuss our results and conclusions in Sec. IV. 
We take h = c = 1 throughout. 

II. EFFECTS ON BIG BANG 
NUCLEOSYNTHESIS 

A time variation in {(p) affects Big Bang nucleosynthe¬ 
sis in several ways. The variation of {(p) changes Gp, SiS 
given in equation CJ, which alters the weak n p inter¬ 
action rates. Changing {(p) also alters the fermion masses 
as 
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Consequently, the electron and the quark masses change. 
The change in the electron mass alters the n ^ p rates 
(see below) and affects the evolution of the electron- 
positron energy density during the epoch of nucleosyn¬ 
thesis. Due to the changes in the u and d quark masses, 
the neutron-proton mass difference changes, and the pion 
mass varies. The latter affects nucleosynthesis through 
the change in the binding energy of the deuteron. We 
have incorporated all of these effects into our calcula¬ 
tions. 

Consider first the neutron-proton mass difference, Q = 
iTin — rup. Following Ichikawa and Kawasaki [^, we take 

Q =-0.76 MeV-b 2.053 MeV^, (3) 

(Wo 

where {4>)o is the present Higgs vacuum expectation 
value. This change in Q alters the ratio of neutron to 
proton abundances in thermal equilibrium: 

^ (4) 

Tip 

The effect on the “^He abundance due to the change in Q 
(alone) is illustrated in Fig.^ An increase in (</>) leads to 
an increase in Q, giving a smaller equilibrium neutron- 
proton ratio, which produces a smaller ^He abundance. 

The weak interactions which interconvert neutrons and 
protons are affected by both the change in Gf and me- 
These interactions are 



FIG. 1: Effects on the '^He abundance (Yp) of a +5% varia¬ 
tion in {(j>) . The solid curve represents Yp as a function of 
the baryon-photon ratio p (7710 = without (0) varia¬ 

tion. Other curves show the isolated effects of changing Gf 
(dot-dash) the neutron-proton mass difference (long dash), 
the electron mass (short dash) and the binding energy of the 
deuteron [with r = 6 in equation Q] (long dash - short dash) 
The dotted curve gives Yp when all the effects are considered 
together. For reference, the observational constraints on Yp 
are shown as two horizontal lines na. 


^ p + e“, 

^ P + Ve, 

^ p + e~ + Pe- 

The total n p rate is 


n + 

n + e'^ 
n 


( 5 ) 
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\{n^ p) = —UnGp / dEe 
^ dme 

{E, + Qf 

1 -b exp[-(ii;e -b Q)/kTi,] 


Ee\Pe\ 

1 -b exp[Ee/kT] 

{Ee - Qyexp{Ee/kT) ) 
1 -b exp[(F;e - Q)/kT^] ) 
( 6 ) 


where the subscripts e and v denote the quantities as¬ 
sociated with the electron and neutrino, respectively. A 
similar expression can be derived for the p —> n rate. 
(Note that in the actual BBN calculation, the weak rates 
are scaled off of the neutron lifetime; we have simply 
rescaled these rates using equation . 

Note that these rates depend on both Gp and rUe, 
both of which are altered by changing {(p) . There is 
a further small change in BBN produced by the effect 
on the expansion rate of the total e+e" energy density, 
which depends on nie- We have included this effect, al¬ 
though it is small. In Fig. 1, we show the effect on the 
^He abundance of changing rrie and Gp separately. Note 
that most of the change from rrie is due to the change in 
the weak rates, as we have noted, rather than from the 


change in the expansion rate. An increase in {(p) results 
in a decrease in Gf, leading to earlier freeze-out of the 
n p reactions, producing more ^He. Similarly, an in¬ 
crease in (p) results in an increase in rrie, decreasing the 
n <-> p reaction rates and also producing more ^He. 

Finally, changing (p) alters the deuteron binding en¬ 
ergy through the change in the pion mass El, 

ml <x ruu + rud <x. (p). (7) 


Early calculations suggested a linear dependence of Bp 
(the deuteron binding energy) on TOtt, and in previous 
discussions of BBN with a time-varying (p) , it was ar¬ 
gued that this effect would be negligible and could be 
ignored Elm (although the effect of changing Bn was 
explored in a slightly different context in Ref. 

Although more recent calculations min have not 
produced a definitive result for the dependence of Bp 
on ttItt, they are in qualitative agreement over our lim¬ 
ited range of interest (i.e., very small changes in ttItt.) In 
this regime. Bp is a decreasing function rriTr (the opposite 
of what was previously assumed El). The calculations 
in both M and Ea have very large uncertainties, but 
within our narrow range of interest, we can safely ap¬ 
proximate the change in Bp as a linear dependence on 
m,r, through 


Bp 

Bdo 


(r-b 1) 


m^ 
— r - 


( 8 ) 
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where the 0 subscript denotes the values of these quanti¬ 
ties at present, and the coefficient r in the fitting formula 
is to be derived from Refs. [3113 We estimate the cen¬ 
tral values for r from these references to be r « 6 and 
r w 10 Although these values for r are quite differ¬ 
ent, we will see that they lead to similar constraints on a 
change in {(j)) . On the other hand, our constraint differs 
sharply from what is derived by ignoring the change in 
Bd- 

The effect on the helium abundance of changing Bd is 
shown in Fig. 1 (for r = 6). It is clear that this effect is 
not negligible. Further, (and unlike the other effects we 
have considered) changing Bd has a significant effect on 
the primordial deuterium abundance. This is most evi¬ 
dent from our graphs of the allowed region for {(j)) (see 
Figs. [21-01 below). An increase in {(p) leads to an in¬ 
crease in m^, resulting in a decrease in Bd- This leads 
to a smaller equilibrium deuterium abundance. Thus, the 
production of ^He begins later, leading to a smaller '^He 
abundance, as shown in Fig. 1. Somewhat paradoxically, 
a decrease in the deuterium binding energy leads to an 
increase in the final deuterium abundance; nucleosynthe¬ 
sis at lower temperatures faces a larger Coulomb barrier, 
leaving more deuterium behind. 

To derive BBN limits on a change in {(j>) , we assume 
two free parameters, the baryon to photon ratio ry, and 
{(j)) / {(j))Q- We vary these parameters within the range 
r] = 10“^° — 10“® and {(p) / (0)g = 0.9 — 1.1. 

To derive our constraints, we consider only the abun¬ 
dances of deuterium and "'^He. Recent observations M 
suggest the limits 


(D/H) = X 10 


-5 


(9) 


on the ratio of deuterium to hydrogen, and 


Yp = 0.238 ±0.005, (10) 

on the primordial mass fraction of "^He. We ignore the 
somewhat more uncertain limits on ^Li in determining 
our constraints, although we have verified that the al¬ 
lowed regions we derive from BBN also produce an ac¬ 
ceptable ^Li abundance. 

Our limits on shown in Figs. El-Q where we have 
translated our constraints on rj into limits on In 

Fig. 121 we show (for reference) the effect of neglecting 
the deuteron binding energy, while Figs. O and 0] cor¬ 
respond to our two choices for r. The solid lines give 
the constraints from the deuterium and ^He abundances, 
with the shaded area indicating the region producing an 
acceptable abundance of both. 

The importance of including the effect of the deuteron 
binding energy is obvious from comparing Fig. [21to Figs. 
Emi In Fig. 121 the regions of allowed deuterium and 
allowed “^He are almost parallel to the axes, indicating 
that the deuterium abundance is nearly independent of 
changes in {(p) ; in this approximation, the abundance of 
deuterium sets the bounds on with the '*He abun¬ 

dance then determining the limits on {(p) . In Figs. rn-TH 



FIG. 2: Regions consistent with the observational data. The 
solid curves represent the constraints from the ^He and deu¬ 
terium abundances, assuming no change in the deuteron bind¬ 
ing energy. The region allowed by BBN is shaded. The dotted 
and the dashed curves represent the 95% and the 68% confi¬ 
dence level regions of the CMB observations, respectively. 





FIG. 3: As Fig. 2, but the BBN limits now assume a change 
in the deuteron binding energy given by equation 0 with 
r = 6. The CMB contours are identical to those in Fig. 2. 


we see a strong dependence of the deuterium abundance 
on {(p) (through the deuteron binding energy), while the 
allowed band for "^He is much narrower. On the other 
hand, Fig. 0and Fig. 0are reasonably similar, allowing 
us to put useful constraints on the time variation in {(p) 
despite the uncertainty in the dependence of Bd on {(p). 

Including the change in the deuteron binding energy, 
we find that the BBN constraint on (cp) variation is ex- 
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FIG. 4: As Fig. 3, with r = 10. The CMB contours are 
identical to those in Figs. 2 and 3. 

tremely tight. Defining A (</>) = (0) — {4>)q, we get 

-0.7 X 10-2 < A(^)/(^)o < 2.0 X 10-2, (11) 

for r = 6, and 

-0.6 X 10-2 < < 1.5 X 10-2, (12) 

for r = 10. 

III. EFFECTS ON THE COSMIC MICROWAVE 
BACKGROUND RADIATION 


FIG. 5: The effect of {(j>} variation on the GMB tempera¬ 
ture anisotrop y p ower spectrum for the ACDM concordance 
model of Ref. IS (Da = 0.66, = 0.12, = 0.2, 

ris = 0.91). The solid line represents the power spectrum 
without {(j)) variation. The dashed and the dot-dashed lines 
correspond to 5% and —5% variations of (0) , respectively. 

where R, (3, Bn, and C represent the recombination co¬ 
efficient, the ionization coefficient, the binding energy of 
the nth hydrogen atomic level, and the Peebles correction 
factor, respectively. The primed quantities stand for the 
modified coefficients following {(j)) variation M- Finally, 
the differential optical depth of the CMB photons is de¬ 
termined by 


Since the weak interaction has no effect on the recom¬ 
bination process, the only effect of changing (^ arises 
through the change in the electron mass, me [ij|. (We 
have neglected any change in the mass of the dark mat¬ 
ter particle. For axions, we expect no change, while the 
change in the mass of a supersymmetric dark matter par¬ 
ticle will be model-dependent). The change in me alters 
both the Thomson scattering cross section and the bind¬ 
ing energy of hydrogen. The Thomson scattering cross- 
section is 


87ra2 2 

while the binding energy of hydrogen is given by 


(13) 



(14) 


The ionization fraction, Xg, is determined by the bal¬ 
ance between photoionization and recombination. The 
evolution of Xg with a variation in ((/>) is given by 


dxg ^ 
dt 


R'm 


2 o'n \ I d3[ B2 

" - /? (1 - a;e)exp (-—— 


( 15 ) 


T = XgUpaT- (16) 

With these modifications, we use CMBFAST to 
generate the theoretical power spectra of temperature 
fluctuations. (See Ref. [ij for a more detailed discus¬ 
sion). Fig. 0 shows the effects of {(p) variation on the 
CMB power spectrum. An increase in {(p) leads to a 
decrease in the Thomson scattering cross section and an 
increase in the hydrogen binding energy. The latter effect 
is the most important, since it shifts the last-scattering 
surface to a smaller comoving length scale, moving the 
curve in Fig. Oto the right. 

Using RAD PACK |^, we perform a analysis for the 
most recent CMB data: VSA, BOOMERANG, DASI, 
MAXIMA and CBI [H H IS El |23- Assuming a 
flat universe (Db + Dcdm + Ua = 1), we set the model 
parameter space 0.5 < ho < 0.8 and 0.3 < Om < 0.4 
(where Dm = Db + ^cdm), and we take the spectral 
index Ug to lie in the range 0.7 < Ug < 1.3. The power 
spectra are obtained by a modified CMBFAST with (p) 
variation. We then marginalize over (p) and . 

The CMB constraints are shown in Figs. I^-TH fidenti- 
cal in all three figures). The dotted and the dashed lines 
represent the 95% and the 68% confidence level regions 
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of the CMB experiments, respectively. In general, the 
CMB calculations provide less stringent constraints than 
BBN; almost the entire region allowed by BBN lies inside 
the 95% confidence region from the CMB. The small area 
allowed by BBN and excluded by the CMB represents a 
tighter bound on but it does not alter the limits 

on {(f)) . 

IV. DISCUSSION 

This work provides new constraints on the possible 
time-variation of ((/>) from both BBN and the CMB. We 
find that the BBN limits are significantly more stringent 
than those that can be derived from the CMB. More 
specifically, we get 

-(0.6-0.7)xl0“^ < A{(j))/{(j)) < (1.5-2.0)xl0“^ (17) 

with the exact numbers depending on the assumptions 
made about the dependence of the deuteron binding en¬ 
ergy on TO^n-. These limits are considerably more stringent 
than those in, e.g., Ref. M , due to better observational 
limits on the primordial element abundances, as well as 
the inclusion of the (previously-neglected) effect of {(j)) 
on the deuteron binding energy. 

In contrast, the limits from the CMB are much less 


stringent. This conclusion is similar to the results ob¬ 
tained for the time-variation in a obtained by Avelino 
et al. [ 5 . (Note also that the effect of a time-variation 
in a and a time variation in {(f>) are almost completely 
degenerate Hi, but BBN breaks this degeneracy). The 
chief advantages of considering the CMB limits are that 
the calculation is much more straightforward (dependent 
only on the change in the electron mass), and that the 
CMB limits should improve sharply in the near future as 
more observational data comes in. Of course, if {(f>) var¬ 
ied between the epoch of BBN (T ^ 10^° — 10® K) and 
the epoch of recombination (T ^ 10® K), then the limits 
we have derived could be applied separately to constrain 
A((/))/((/)) at each epoch. 
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